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ABSTRACT. We investigate the properties of the collision operator Q associated to the linear Boltz- 
mann equation for dissipative hard-spheres arising in granular gas dynamics. We establish that, as 
in the case of non-dissipative interactions, the gain collision operator is an integral operator whose 
Qi^ ' kernel is made explicit. One deduces from this result a complete picture of the spectrum of Q in an 

Hilbert space setting, generalizing results from T. Carleman (6) to granular gases. In the same way, 
we obtain from this integral representation of Q + that the semigroup in L 1 (R 3 X R 3 , dx <8> dv) asso- 
ciated to the linear Boltzmann equation for dissipative hard spheres is honest generalizing known 
results from (T|. 

KEYWORDS. Granular gas dynamics, linear Boltzmann equation, detailed balance law, spectral 
theory, Co-semigroup. 



J> ! 1. Introduction 

m 

We deal in this paper with the linear Boltzmann equation for dissipative interactions modeling 
the evolution of a granular gas, undergoing inelastic collisions with its underlying medium. Actu- 
ally, we shall see in the sequel that there is no contrast between the scattering theory of granular 
C — > ■ gases and that of classical (elastic) gases. This may seem quite surprising if one has in mind the 

fundamental differences that may be emphasized between the nonlinear kinetic theory of granular 
_C 1 gases and that of classical gases, as briefly recalled in the next lines. 

03 ' 

1.1. Granular gas dynamics: linear and nonlinear models. Let us begin by recalling the gen- 
eral features of the kinetic description of granular gas dynamics that can be recovered from the 
monograph [4] or the more mathematically oriented survey (23J. If f(x,v,t) denotes the distribu- 
^ ■ tion function of granular particles with position x € R 3 and velocity v € R 3 at time t > 0, then the 

evolution of f(x, v, t) is governed by the following generalization of Boltzmann equation 

d t f(x, v,t) + v V x f(x, v, t) = C(f)(x, v, t), (1.1) 

with initial condition f(x, v, 0) = fo(x, v) 6 L X (R 3 X R 3 , dx ® dv), where the right-hand side C(f) 
models the collision phenomena and depends on the phenomenon we describe. 

In the nonlinear description, the collision operator C(f) -: S[f,f] is a quadratic operator 
modeling the binary collision phenomena between self-interacting particles. For hard-spheres 
interactions, it reads 



mfl(v)= f \q-n\ 

Jr 3 xS 2 



1 

— f(x, v+, t) f(x, w+, t) - f(x, v, t) f(x, w, t) 



dwdn, 
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where q is the relative velocity, q — v — w. The microscopic velocities (v+,w+) are the pre- 
collisional velocities of the so-called inverse collision, which results in (v, w) as post-collisional 
velocities. The main peculiarity of the kinetic description of granular gas is the inelastic character 
of the microscopic collision mechanism which induces that part of the total kinetic energy is 
dissipated. This energy dissipation might be due to the roughness of the surface or just to a non- 
perfect restitution and is measured through a restitution coefficient < e < 1 (which we assume 
here to be constant for simplicity, see Remark [L2l ). As a consequence, the collision phenomenon 
is a non microreversible process. Generally, we assume that the energy dissipation does not affect 
the conservation of momentum. Therefore, in the homogeneous setting, i.e. when fo(x, v) - fo(v) 
is independent of the position, the number density of the gas is constant while the bulk velocity is 
conserved. However, the temperature of the gas 



continuously decreases (cooling of granular gas). As a consequence, the stationary state of the 
inelastic collision operator S is a given Dirac mass. However, the homogeneous Boltzmann equa- 
tion for granular gases exhibits self-similar solution (homogeneous cooling state) J5J[T7l. Note 
the important contrast with the classical kinetic theory, i.e. when e = 1, for which it is well-known 
that the steady state of the collision operator is a Maxwellian distribution. 

The linear Boltzmann equation for dissipative interactions concerns dilute particles (test 
particles with negligible mutual interactions) immersed in a fluid at thermal equilibrium lfl4l IT6l 
l22l . The total kinetic energy is dissipated when the dilute particles collide with particles of the host 
fluid. Such physical models are well-suited to the study of the dynamics of a mixture of impurities 
in a gas (9j HI for which the background is in thermodynamic equilibrium and that the polluting 
particles are sufficiently few. We refer the reader to [11] and the survey ifTUl for more details on 
the theory of granular gaseous mixtures. Assuming the fluid at thermal equilibrium and neglecting 
the mutual interactions of both the test and dilute particles, the collision operator C(f) = Q(f) is 
a linear scattering operator given by 



where Ai\ stands for the distribution function of the host fluid. Note that in such a scattering 
model, the microscopic masses of the dilute particles m and that of the host particles m~\ can be 
different. We will assume throughout this paper that the distribution function of the host fluid is a 
given normalized Maxwellian function: 



where U\ € R 3 is the given bulk velocity and $i > is the given effective temperature of the host 
fluid. It can be shown in this case that the number density of the dilute gas is the unique conserved 
macroscopic quantity (as in the elastic case). The temperature is still not conserved but it remains 
bounded away from zero, which prevents the solution to the linear Boltzmann equation to converge 
towards a Dirac mass. This strongly contrasts to the nonlinear description and suggests that the 




Q(f) = S[f,Mi]= f \q-n\\\f(x / v ir/ t)M 1 (w ir )-f(x / v / t)M 1 (w)]dwdn (1.2) 
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linear scattering model associated with granular gases does not contrast too much with the one 
associated with classical gases. 

The first mathematical result in this direction is the following one according to which, as in the 
classical case, the unique steady state of Q remains Gaussian. The fact that the linear Boltzmann 
equation still possesses a stationary Maxwellian velocity distribution was first obtained in [ 16 ] and 
we refer to Ifl4l for a complete proof (existence and uniqueness) for hard-spheres model (see also 
E2l for a version of this result for Maxwell molecules) : 

Theorem 1.1. The Maxwellian velocity distribution: 

. „, ^ I m \ 3 / 2 f m(v-ui) 2 ) o 

# (1 + e)m 

with 9 = ; - $1 is the unique equilibrium state ofQ with unit mass. 

2m + (3 + e)m 1 H H J ^ 

Remark 1.2. Note that, if one does not assume the restitution coefficient e to be constant (see 
(4J for the general expression of non-constant restitution coefficient e — e(q) in the case, e.g., of 
visco-elastic spheres) then the nature of the equilibrium state of Q is still an open question: it is 
not known whether such a steady state is a Maxwellian or not. Consequently, it is still not clear 
that linear inelastic scattering models behave like elastic ones. For this reason, we shall restrict 
here our study of the linear Boltzmann equation to a constant restitution coefficient. We also 
point out that, if the distribution function of the host fluid M\ is not of gaussian type, the explicit 
expression of the equilibrium state of Q is an open question to our knowledge. 

The existence and uniqueness of such an equilibrium state allows to establish a linear version 
of the famous H-Theorem. Precisely, for any convex C 1 -function O : R + — > R, one can define 
the associated entropy functional as 

Theorem 1.3 (H-Theorem lfT4ll20l ). Let fo(v) be a space homogeneous distribution function with 
unit mass and finite entropy, i.e. Hq,(fo\Af) < 00. Then, 

jH„(f(t)\M) < (t > 0), (1.4) 

where fit) stands for the (unique) solution to (11.11 ) in L 1 (R 3 ,dv). 

Note that such a result is valid for any scattering operator with positive kernel and positive 
equilibrium ||20l . As an important consequence, it can be shown by suitable compactness argu- 
ments that any solution to the Boltzmann equation dl.lb (with unit mass) converges towards the 
Maxwellian equilibrium At. Note also that, for the nonlinear Boltzmann equation for dissipative 
interactions, the temperature is a trivial Lyapunov functional leading to the convergence of any so- 
lution towards a delta mass. However, the construction of a Lyapunov functional in the self-similar 
variables allowing relaxation towards the homogeneous cooling state is still an open question (see, 
e.g. ifTTl for related problems). 



4 



L. ARLOTTI & B. LODS 



To summarize, the steady state of the linear collision operator for dissipative interactions has 
the same nature (a Maxwellian distribution) as the one corresponding to non-dissipative interac- 
tions. Moreover, as in the classical case, by virtue of the H-Theorem, such a steady state attracts 
any solution to the space homogeneous Boltzmann equation dl.ll) . This seems to indicate that 
most of the properties of the linear Boltzmann equation for elastic interactions remain valid for 
inelastic scattering models. It is the main subject of this paper to make precise and confirm such 
an indication and the key ingredient will be the derivation of an integral representation of the gain 
part of the collision operator. 

1.2. Main results. The main concern of our paper is the derivation of a suitable representation of 
the gain part of the collision operator Q as an integral operator with explicit kernel. Precisely, the 
linear collision operator Q can be split into Q(f) = Q + (f) — Q~(f), where the gain part is 

Q + (f)(v) = e- 2 f Iq-nlfivJMtiwJdwdn 

Jr 3 xs 2 

while 

Qif)( v ) = I l<? ' n\f(v)M\{w)dwdn = o(v)f(v) 

J«3xS 2 

where the collision frequency a(v) is given by a{v) — I*? ' n\M\(w)dwdn. It is well-known 

that, for non-dissipative interactions, i.e. when e = 1, the gain part Q + can be written as an 
integral operator with explicit kernel (H [T5J (see also lfl2l 13 for similar results for the linearized 
Boltzmann equation). We prove that such a representation is still valid in the dissipative case: 

Theorem 1.4. Iff > is such that o(v)f(v) £ L^IR 3 ,^), then 

Q + (f)(v)= f k{v,v')f{v')dv> 

JR3 

where the integral kernel k{v,v') can be made explicit (see (12.21 ) ). 

Actually, most important is the fact that the integral kernel k(v,v') turns out to be very similar 
to that obtained in the classical case (see for instance |[T5l l6l). the only changes standing in some 
explicit numerical constants. Moreover, as we shall see, the kernel k(v,v') and the Maxwellian 
distribution M satisfy the following detailed balance law: 

k(v, v')M{v') = k{v', v)M{v), v, v' e R 3 , 

that allows us to recover Theorem ll.ll in a direct way. Recall that, in iTfll . the Gaussian nature of 
the steady state of Q was obtained by replacing Q by its grazing collision limit. 

We derive from these two results some important consequences on the linear Boltzmann equa- 
tion (11.11 ) with C = Q. The applications are dealing with the space dependent version of (11.11 ) as 
well as with the space homogeneous version of it. The first one concerns the spectral properties of 
the Boltzmann collision operator in its natural Hilbert space setting. 
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1.3. Spectral properties of the Boltzmann operator in L 2 (At _1 ). Applying the above H-Theo- 
rem |1.3| with the quadratic convex function O(x) = (x — l) 2 , one sees that a natural function space 
for the study of the homogeneous linear Boltzmann equation is the weighted space L 2 (yM _1 ). 
Now, from Theorem 11.41 it is possible to prove that the gain collision operator Q + is compact in 
L 2 (A1 _1 ). This compactness result has important consequences on the structure of the spectrum 
of Q as an operator in L 2 (A1 _1 ). Precisely, from Weyl's Theorem, the spectrum of Q in this space 
is given by the (essential) range of the collision frequency a(-) and of isolated eigenvalues with 
finite algebraic multiplicities. Since A = is a simple eigenvalue of Q (its associated null space is 
spanned by Af), this leads to the existence of a positive spectral gap. In turns, one proves that any 
solution to the space-homogoneous linear Boltzmann equation (11.11 ) converges at an exponential 
rate towards the equilibrium. These spectral results are technical generalizations of some of the 
fundamental results ofT. Carleman [6], but are new in the context of granular gas dynamics. 

1 .4. Honest solutions for hard-spheres model. It is easily seen that, for any nonnegative /, 

f Q + (f)(v)dv= f o{v)f{v)dv, (1.5) 

JR3 JR 3 

i.e. the collision operator Q is conservative. Then, formally, any nonnegative solution f(x, v, t) to 
(11.11 ) (with C — Q) should satisfy the following mass conservation equation: 

I f(x,v,t)dxdv = I f(x,v,0)dxdv, Vf > 0. (1.6) 

JR 3 xR 3 Jr 3 xR 3 
It is the main concern of Section |4] to prove that such a formal mass conservation property holds 
true for any nonnegative initial datum f(x, v, 0) € L X (R 3 X R 3 ). As well documented in the mono- 
graph [3], this is strongly related to the honesty of the Co-semigroup governing Eq. (ll.ll ). More 
precisely, if we denote by To the streaming operator: 

@(T ) = {feX,v- V x / € XL T / = -v ■ V x f, 

it is not difficult to see that there exists some extension G of To + Q that generates a Co-semigroup 
of contractions (Z(t)) t ^o in X = L 1 (R 3 X R 3 ). According to the so-called "sub-stochastic pertur- 
bation" theory, developed in Hl|3]|25]], it can be proved that 

I Z(t)f(x,v)dxdv = I f(x,v)dxdv, V/eX,/>0 

Jr 3 xR 3 ' Jr 3 xR 3 

if and only if G is the closure of the full transport operator: G — Tq + Q. We show in Section 4 that 
the latter holds. To do so, we shall use the integral representation (Theorem 1 1.41) in order to apply 
some of the results of (2 (see also [3, Chapter 10]) dealing with the classical linear Boltzmann 
equation. 

1.5. Organization of the paper. We derive in Section 2 the integral representation of Q + (The- 
orem 12.11) as well as some of its immediate consequences concerning the explicit expression of 
the collision frequency. We also recover Theorem 11.11 throu gh a detailed balance law. Section 3 is 
devoted to the study of the collision operator Q in the narrow space L 1 (M~ 1 (v)dv) and its spectral 
consequences. In Section 4 we apply the results of Section 2 as well as some known facts about 
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the classical linear Boltzmann equation H] O to the honesty of the solutions to the Boltzmann 
equation for dissipative hard-spheres. 



2. Integral representation of the gain operator 
Let us consider the gain operator for dissipative hard-spheres: 

Q + (/)(y) = e- 2 f \q-n\f{v^)Mi{w^)dwdn 

jR3xS 2 

and let o(v) be the corresponding collision frequency: 

a(v) -I \q ■ n\Mi{w)dwdn, »6R 3 . 

Jr 3 xS 2 

Recall that M\ is a Maxwellian distribution function with bulk velocity U\ and effective tem- 
perature $i. We recall here the general microscopic description of the pre-collisional velocities 
(v+, w+) which result in (v, w) after collision. For a constant restitution coefficient < e < 1, one 

has aiia 

1-6 

= p- 2a 1 [q- n]n, 
1-B 

w+ = w + 2(1 - oc) l _ 1 p [q ■ n]n; 
where q — v - w, a is the mass ratio and $ denotes the inelasticity parameter 

oc = , £ = — — . 

m + m-[ 2 

We show in this section that, as it occurs for the classical Boltzmann equation, Q + turns out 
to be an integral operator with explicit kernel. The proof of such a result is based on well-known 
tools from the linear elastic scattering theory (6j [12l EES while, in the dissipative case, similar 
calculations have been performed to derive a Carleman representation of the nonlinear Boltzmann 
operator in iTTTl . 

Theorem 2.1 (Integral representation of Q + ). For any f e L^R 3 x R 3 , dx <g> o{v)dv), 

Q + f(x,v) = -t ri [ f{x,v')k{v,v')dv' , (2.1) 
2e 2 y 2 Jr3 



where 



KV ' V) = \2rth) lV ~ Vl exp -— ^ ( 1 + M) |,-,| + — j (2.2) 



2a(l - JS) - 1 1 — jS 

with n = — — — > and y - a-. 



a(l-/3) 1 1-2)3 
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Proof. The local (in x) nature of Q + is obvious and we can restrict ourselves to prove the result 
for a function / € L l {E? ,o{v)dv) that does not depend on x. Set y - oc^J^ and y = (1 - «)t~^ 
so that 

v+ = v- 2y[q ■ n]n and -w + 2y[q- n]n. 
The following formula, for smooth (p: 

f (q • n)(p ((q • n)n) dn = ^ f <p\- — \do = \ f 5(2x ■ q + x 2 )(p(-x/2)dx. 
Js + 4 Js 2 \ 2 / 2 Jr 3 

applied to 

cp(x) - f(v- 2yx) Mi (w + 2yx) 

yields 

<TM = e-* f 5(2,-, + ^)/(, + rx)At l(l ,-rx)^. 

jR 3 xR 3 

The change of variables x i-> v' - v + yx leads to 

Q + f(v) = e~ 2 y- 3 \ 5(2y- 1 (v' -v)-q + y~ 2 \v' - v\ 2 )f{v')Mi{w - ~{v' - v))dwdv' . 
Jr 3 xr 3 ' 7 

y 

Now, keeping v and v fixed, we perform the change of variables w h-> w - w (v' - v), which 

leads to 

Q + f(v) = e~ 2 y- 3 J 5\2y- 1 (v' - v) ■ [v - w' - Y -{v' - v)\ + y~ 2 \v' - v\ 2 \ x 

JR3XR3 \ Y ' I 

xf(v')Mi(w')dw'dv'. 

Writing w' - v + X\n + V2 with A\ = (w' - v) • n e R, n = (v' - v)/\v' - v\ and V2 ■ n - 0, we 
get, noting that dw' - dV^Ai , 

Q + f(v) = e" 2 7" 3 f f(v')dv' [ d\ x f M 1 (v + V 2 + A 1 n)dV 2 x 
Jr 3 Jr JV r n=0 

x 5 (y"V - v\ 2 - 2yy~ 2 \v' - v\ 2 - 2y~ 1 A 1 \v' - v\) . 
Thanks to the change of variables A\ i-> 2y~ 1 \v' - v\A\, one can evaluate the Dirac mass as 

f 6 (y"V - v\ 2 - 2yy- 2 \v' - v\ 2 - 2y~ 1 A 1 \v' - v\) Mi(v + V 2 + A x n)dA x 
Jr 

Mi [v + V 2 + -^—^-(v' - v) 



2\v'-v\ l \ 2y 
where we used that n = (1/ — v)/\v' - v\. Consequently, 



Q + f{v) = ^^k{v,v')f{v')dv' 
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where 



k(v,v') = —^— f M 1 {v + V 2 + ^—^(v'-v))dV 2 . 

\v ~ v '\ Jv 2 iv'-v)=0 \ 2 Y ) 

It remains now to explicit k(v,v'). We will use the approach of lfl5"1 . Let us assume v,v' to be 
fixed. Let P be the hyperplan orthogonal to (v' - v). For any V 2 e P, set 

l-2y , 
z = v + ^ ( v -v) + V 2 -u 1 

so that 

k(v,v') = (^] \v-vT 1 f exp{- Ql z 2 }dV 2 . 

mi v + v' 1 — 2y 
where Q\ = Denoting for simplicity u - — U\ and /j = , one has 



z 2 



= l« + v 2 \ 2 + ^^-b - *f + ^(1* - "il 2 - W - i/il 2 ) 

where we used the fact that V 2 is orthogonal to (v' - v). Splitting u as 

U - Uq + Mj_ 

where Uq is parallel to P — v' while w ± is orthogonal to V — v' (i.e. u± 6 P), we see that 

\\v - Ul \ 2 - \v' - u^ 2 ] 

\u + V 2 \ 2 = \u \ 2 + \u ± + V 2 \ 2 and |m | 2 = -= 

so that 



2 



fcfo t/) = \v- v'r 1 (^) 3/2 exp (- ei \u ± + y 2 | 2 ) rfy 2 



ex Pi-j 



\\v - U\\ - \V -U\' 1 

(1 + y) 2 \v - v'\ 2 + 2(1 + n)(\v - «i| 2 - - wil 2 ) + — 



Finally, since u± e P, 

I exp(-pi|w_L + ^ 2 | 2 W 2 = I exp(-^x 2 )dx = — , 
Jp v 7 Jr 2 01 

one obtains the desired expression for k(v, v'). □ 

The very important fact to be noticed out is that the expression of k(v,v') is very similar to 
that one obtains in the elastic case [15], the only change being the expression of the constant /i. 
In particular, in the elastic case e = 1, we recover the expression of the kernel obtained in [6] for 
particles of same mass (i.e. m = m{) and in [ 15 ] for particles with different masses. 
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Another fundamental property of the kernel k(v,v') is that it allows us to recover the steady 
state of Q through some microscopic detailed balance law. Precisely, 

Theorem 2.2. With the notations of the Theorem \2.1\ the following detailed balance law: 

k(v, v') exp + ~ "i) 2 } = W, v) exp (~d + ~ "i) 2 } , (2-3) 

holds for any v, v' 6 IR 3 . As a consequence, the Maxwellian velocity distribution: 

/ m v3/2 f m (p-Ui) 2 \ 3 

, (l-a)(l-/3) 

with $ - — — — * s ^ e u ™ aue equilibrium state of Q with unit mass. 

Proof. According to Eq. (12.21 ), it is easily seen that 

k{v',v) = k(v,v')exp{^(l + ^) - Ml | 2 - \v' -«i| 2 )}, v,xt e R 3 

which is nothing but (12- 3b - Now, writing ^-(1 + fi) = , straightforward calculations lead to the 

desired expression for the equilibrium temperature 9*. The fact that M is an equilibrium solution 
with unit mass follows then from the fact that 

Q(M)(v) = f k(v, v')M(v')dv' - o{v)M{v) = f [k(v, v')M(v') - k(v', v)M{v)} dv' 

J«3 JR3 



and from the detailed balance law (12.3b - To prove that the steady state is unique, we adopt the 
stategy of |[2T1 Theorem 1]. Precisely, consider the equation 

o(v)f(v) = Q + f(v), V^elR 3 (2.4) 

which admits at least the solution / = At. Since o(v) does not vanish, any solution / to (12.41 ) is 
such that 

/(l,) = ^) Q+(/)(y) ' Vy£R3 ' 
Since Q + is an integral operator with nonnegative kernel, it is clear that o(v)\f(v)\ < Q + (|/|)(^) 
for any »eR 3 . Now, from the positivity of both a and Q + , one sees that the conservation of mass 
(11.51 ) reads: 

\\°f\\x = \\°\f\\\x = \\Q + {\f\)\\x. 
This shows that, actually, \Q + (f)(v)\ = a{v)\f(v)\ - Q + (\f\)(v) for any v e R 3 . Again, since Q + is 
a positive operator, one obtains that 

f=m 

Now, assume that (12.41 ) admits two solutions f\, fi with f 3 fi(v)dv = J R3 f2(v)dv - 1. Then, 
fx - fi is again a solution to ( 12.41 ) so that, f\— f% = ±|/i - fi\- Thus, 

± r i/i(y)-/ 2 (*= r /i^do- r ^^=0 

JR3 JR3 JR3 

and the uniqueness follows. □ 
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The above result allows to derive the explicit expression of the collision frequency a(v): 
Corollary 2.3. The collision frequency o{v)for dissipative hard-spheres interactions is given by 

2n I m x (A&t ( m x 2 \ 

+ ( 2 \ v - Ul \ + -^-\ r~ U \ xv (-^A d X (2.5) 

Consequently, there exist positive constants vq, V\ such that 

v (l + \v - «i|) < a{v) < vi(l + fa - Mil), Vw £ R 3 . 

Proo/ Set C = ^ ■ Noting that a(v) - J K3 k(v' , v)dv' for any v € R 3 , one has, with the 
change of variable z = v' - v, in a polar coordinate system in which v lies on the third axis 

,\ nC f m l (n , m i b-wil 2 -|2: + y-«il 2 \ 2 l , ,_i , 
a(,) = C j R3 exp|-— + ^ J |N dz 

= 2nC I dp I gexpl-— — ((2 + ^l)q + 2\v - Ui\ cos cp) z \ sin q) dep. 
Jo Jo I 8#i ' " J 

The computation of this last integral leads to the desired expression for o(v). The estimates are 

then straightforward iPTSl . □ 

3. Application to the Boltzmann operator in L 2 (M~ l ). 
We investigate in this section the properties of the Boltzmann operator Q in the weighted space 

<H = L 2 (R 3 ;M~Hv)dv). 

We shall denote by the inner product in < H. The introduction of such an Hilbert space setting 

is motivated by the application of the H-Theorem 11.31 with the convex function O(x) — (x — l) 2 . 
In this case, one sees that, if /o > is a space homogeneous initial distribution such that 

f f (v)dv = l, f ^{vfM^WvKoc, 

JR3 Jr3 
then any solution f(t, v) to the space homogeneous equation 

dtffrv) = Q{f){t,v), f(0,v) = f (v) € <H, (3.1) 
satisfies the following estimate: 



t > 0. 



In other words, the mapping t i — > •) - At||^ is nonincreasing. For these reasons, the study 
of the properties of the collision operator Q in *H is of particular relevance for the asymptotic 
behavior of the solution 

d t f(t, v) = Q(f)(t, v), /(0, v) 6 7/. (3.2) 
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The material of this section borrows some techniques already employed by T. Carleman |6] in the 
study of non-dissipative gas dynamics (see also, e.g. Q or |[T2l for similar results in the context 
of the linearized Boltzmann equation). Let _£ be the realization of the operator Q in < H, i.e. 



J \f(v)\ 2 o(v)M-Hv)dv< ooj. 



and, for any / € &{■£.), -Lf{v) = Q(f)(p) is given by (11.2I ). As previously, one can use the 
following splitting of X as a gain operator and a loss (multiplication) operator, X = -C + - -C~ with 

£ + (/)(z;) = f k{v,v')f{v')dv' and £~{f) = o(v)f(v), f e 9){£). 

We shall show, as in the classical case, that X + is actually a bounded operator in f H. Precisely, let 
^ define the natural bijection operator from L 2 (R 3 , dv) to 7Y: 

J : L 2 (R 3 , dv) — > H 

It is clear that is a bounded bijective operator whose inverse is given by 
f- l g{v) = M- 1/2 (v)g(v) € L 2 (R 3 , cfc), Vg € <H. 

Now, let us define 

G(v, v') = M~ 1/2 (v)k(v, v')M 1/2 {v'), v, v' e R 3 , 

i.e. 

r< >\ ( m i \ 1/2 | 'i-i / m WnWi /|2 (jp - u x | 2 - \v' - urf) 2 ]] 

(3-3) 

From the detailed balance law (12.31 ). one easily checks that G(v, v') = G(v', v) for any v, v' e R 3 . 
Therefore, defining Q as the integral operator in L 2 (R 3 , dv) with kernel G(v, v'), i.e. 



em= f G{v,v')f{v')dv' , 

Jr 3 



one can prove the following: 



Proposition 3.1. Q is a bounded symmetric operator in L 2 (R 3 , dv) and £ + = J? Q J? . Conse- 
quently, X + is a bounded symmetric operator in ^H. 

Proof. It is clear that Q is symmetric since G(v, v') = G(v', v). Now, to prove the boundedness of 
Q, one adopts a strategy already used in the non-dissipative case by T. Carleman [6, p. 75] and 
shows easily that 

C := sup I G(v,v')dv' < oo. 

veK? Jk 3 
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Since G(-, •) is symmetric, one also has sup^,^ J R3 G{v,v')dv — C < oo. Denoting by (•, •) the 
usual inner product of L 2 (R 3 ,dv), one deduces from Cauchy-Schwarz identity, 

Iff \f(v)\ 2 dv+ f Wfdv'), Vf,geL 2 (R\dv), 
1 \Jr3 Jr 3 / 

which leads to the boundedness of Q. Now, since G{v,v') = M~ 1 ^ 2 (v)k(v,x/)Ai 1 ^ z (v / ) for any 
v, v' e R 3 , one gets easily that X + = an d the conclusion follows. □ 

In Proposition 13.11 we proved that the gain operator X + is bounded in < H, i.e. X + € SSi^H). 
Actually, we have much better and it is possible, as in the non-dissipative case, to prove that £ + 
is a compact operator in < H. Precisely, the following lemma is a direct consequence of Theorem 
12. II and similar calculations valid for the non-dissipative case [6, p. 70-75]. However, we give a 
detailed proof of it since the known similar results by T. Carleman are all dealing with the case 
m = m\ and e = 1. It has to be checked that taking account the parameters m rtl\ and e < 1 
does not lead to supplementary difficulty (see Remark [3T8l where the role of e ^ 1 does not allow 
to adapt mutatis mutandis a result valid in the elastic case). 

Lemma 3.2. For any < p < 3 and any q > 0, there exists C(p, q) > such that 

Proof. The proof is a technical generalization of a similar result due to T. Carleman [6] in the 
classical case (i.e. when m — m\ and e = 1). Let us fix < p < 3 and q > and set 

Then, one sees easily that 

/ m \P/2 r n r°° expj-^((l + ^)V + (? + 2|^-Mi|cos(p) 2 )} 

I(p)= 2n(^ r ) smcpdcpi q 2 ~v 1 8 ^ ^ —do. 

\2n$ x ) Jo Jo (l + ^q 2 + \v- Mil 2 + 2q\v - ui|cos<p] 

Note that, since < p < 3, 

sup I(v) < 471 {^-f 2 [ Q 2 ~ p exp {-^(1 + tfg 2 ) d 6 < oo. (3.4) 
veE? \2n9\l Jo I 8-Si J 

Performing the change of variable x = - U\\ + 2cos(p, y = q/\v - U\\, one has (x, y) € Q 
where 

Q = {(x, y) € K 2 ; y > 0, |x - y|< 2} 

and 

12 ^ dxdy 



. . , /mi \P /2 . l3 _ p f f m x ip\v-u^- , 2 2 2 v 



yP~ 2 (l + \v~ui\ ^1 + xy) 



v 
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We split Q into Q = Qi U Q2 where Qi is the half-ellipse 

Qi = {{x, y) € R 2 : y > 0, (1 + \i) 2 y 2 + x 2 < 1/4} while Q 2 = Q \ Qi. 

Note that, since 1 + y. > 1, one has Qi C Q. One defines correspondingly I\(v) and l2(v) as the 
above integral over Qi and Q2 respectively. One notes first that, if (x, y) € Qi then xy > - g^ + ^ 
so that 

1 - 8 ^+„) , < a < 1. Letting R = (f^-) and setting f = _R|u - U\\x, u = 

R(l + f0|t>-«lly, it is easy to check that -R\v~U\\l2 < t < R|w-Mi|/2, while < u < R|w-Ui|/2, 
so that 

U '"l2nsW b-«il(l+*-ml)' J K J„ ul- 2 
Thus, there exists a constant Ci(p, (?) > such that 

^ < 1 iTTTT iv?' Vy 6 R3 - (3 " 5) 

|u - «i|(l + - Wl|) 9 

Let us now deal with hip). Arguing as above, 

IM - (^f -I" - ".I 3 "" £ ex P ((1 + ,)V + x 

x — dxdy. 

yP" 2 (l + |o - ui\ yjl + xy] 

Clearly, since (1 + /i) 2 y 2 + x 2 > 1/4 for any (x, y) € Q2, then 

,1 dxdy 



f-2 



I—) ?exp( — exp( _( 1+(j)y j_j 



Hence, there is some constant C^{p, q) such that 

Hv) < C 2 (p, 9 ) exp |_ R2 l p -"il 2 j , p e R 3. (3 . 6) 
Combining (13.51 ) and (13.61 ). one sees that 

M«i ,nff ^ + C ^^)exp(-g^), 

\v - Ui\(l + a\v - \ 8 / 

According to (I3.4I ). lim sup| B _ M < °°/ from which we get the conclusion. □ 
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Remark 3.3. Note that the above Lemma can be extended to more general collision kernels ( in- 
cluding long-range interactions) following the lines of the recent results [18 ] dealing with the 
elastic case. 

From the above Lemma, one has the following compactness result: 

Proposition 3.4. Q is compact in L 2 (M. 3 ,dv). Consequently, £ + is a compact operator in < H. 

Proof. Applying arguments already used in (6), the above Lemma implies that the third iterate of 
Q is an Hilbert-Schmidt operator in L 2 (M?,dv), i.e. the kernel of is square summable over 
E 3 X R 3 . The compactness of Q follows then from standard arguments and that of X + is deduced 
from the identity X + = J? Q J*~ x (see Proposition [3J}. □ 

The following, which generalizes a known result from classical kinetic theory, proves that X is 
a negative symmetric operator in IH: 

Proposition 3.5. The operator {£, &{£)) is a negative self— adjoint operator ofH. Precisely, 

n 2 



<£/,/>* = 4 f k{v,v')M{v')\M- l {v)f{v)-f{v')M- l {v') 



dvdv' «S 



for any f € &(£). 

Proof. The fact that (X, is self-adjoint is a direct consequence of Proposition 13. ll since X~ 

is clearly symmetric. Now, it is a classical feature, from the detailed balance law (12.3I ). that 

(£f,f) n = f k{vM)M{v')\M- l {v')f{v')- f^M-Hv^f^M^^dvdv' . 

Exchanging v and v' and using again the detailed balance law (12.31) . one sees that 

{£f r f)^ = [ k(v,v')M(v')\M-\v)f(v) - f(v')M-Hv')]f(v')M-\v')dvdv' 

JR3XR3 L -I 

so that, taking the mean of the two quantities, 



{Lf,f) n = -\ f k{v,v')M{v')\M-\v)f{v)- f{v')M-\v') 

L JR3XR3 L 



dvdv' < 



which ends the proof. □ 

Remark 3.6. From the above result, the spectrum S(X) of £ lies in R_, i.e. S(X) c (-00,0]. It 
is clear that A = lies in S(X). Precisely is a simple eigenvalue of £ since M is the unique (up 
to a multiplication factor) steady state of £. 

Combining the above results with Proposition 13. ll leads to a precise description of the spectrum 
of X: 

Theorem 3.7. The spectrum of £ (as an operator in 1H) consists of the spectrum of —£~ and of, 

at most, eigenvalues of finite multiplicities. Precisely, setting vq — infj, sK 3 a{v) > 0, 

S(X) = {A £ R; A < -v } U(A„;«6J| 

where I C IN and (A n ) n is a decreasing sequence of real eigenvalues of £ with finite algebraic 
multiplicities: Aq = > Ai > A2 • • • > A M > . . . , which unique possible cluster point is —vq. 
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FIGURE 1 . Spectrum of the collision operator Q in *H. 

Proof. From Proposition 13.41 X is nothing but a compact perturbation of the loss operator —£r. 
Hence, Weyl's Theorem asserts that S(X) \ S(— -£~) consists of, at most, eigenvalues of finite 
algebraic multiplicities which unique possible cluster point is sup{A, A 6 S(— JL~)}. In particular, 
up to a rearrangement, one can write S(-£) \ ®(— £~) — {A n/ n € 1} with Ao > A\ > Ki ■ ■ ■ > 

A n > We already saw that Ao = since M is a steady state of Q and Ai € 'H. Now, since 

—£T is a multiplication operator by the collision frequency -a(-), its spectrum S(-X~) is given 
by the essential range R ess (-a(-)) of the collision frequency From Corollary 12.51 one sees without 
difficulty that 

Ress(-o(-)) = (~°°/ -vo] 
where v = inf oe]R 3 cr(p) = lim| o _ Ml |^ o(v) = is positive. □ 

Remark 3.8. We conjecture that, as it is the case for elastic interactions 1131 . sef of eigenvalues 
lying in (—vq, 0) infinite. However, the technical generalization of the proof of 11131 appears to 
be non trivial because of the non zero parameter /i. We thank anyway an anonymous referee for 
having pointed to us the reference lfl3l . 

The above result provides a complete picture of the spectrum of Q as an operator in 'H (see 
Fig. 1) and shows, in particular, the existence of a positive spectral gap |Ai| of £,. Note that such 
a result, combined with Proposition 13.51 has important consequence on the entropy production, 
since it can be shown in an easy way that the H-Theorem reads as 

j t \\f{t)-Mf H = (m\f(t))<H- 

Consequently, the Dirichlet form B(f) = (_£/ : ,f)<H pl a Y s the role of entropy-dissipation functional 
and the existence of a spectral gap |Ai| is equivalent to the following coercivity estimate for £>(/): 

B(/)>-|Ailll/"||^ V/±span(M). 
One deduces easily the following corollary on the exponential trend towards equilibrium: 

Corollary 3.9. Let fo(v) € *H and let f(t,v) be the unique solution to the linear homogeneous 
Boltzmann equation (13 -2b - Then, there is some constant C > such that 

\\f(t,-)-M\\„<Cexp(-\Ai\t)\\f -M\\„, for any t > 0, 

where < |Ai| < vq is provided by Theorem \3.7\ 
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We refer the reader to UM for details on the matter, and in particular, for an explicit estimate 
of the spectral gap \A\\. 

4. Application to the honest solutions of the Boltzmann equation 

4. 1 . Conservative solutions. We are interested in this section in applying the result of the previ- 
ous section to prove the existence of honest solutions to the linear Boltzmann equation for dissi- 
pative hard-spheres 

d t f(x,v,t) + v ■ V x f{x,v,t) = Q{f){x,v,t), (4.1) 

with initial condition 

f(x, v, 0) = f (x, v) € L J (R 3 x R 3 , dx <8> dv), 

where the collision operator Q is given by Eq. dl.2| ). Recall that the streaming operator To is 
defined by 

9(T ) = {feX,v-V x feX}, T f = -v-V x f 

where X = L^R 3 X R 3 , dx ® dv). One can define then the multiplication operator E by 

9(L) = {feX,ofeX}, E f{x, v) = -a(v) f(x, v) 

where, as in the previous Section, a(v) is the collision frequency corresponding to dissipative hard 
spheres interactions and given by Eq. (I2.5I ). The following generation result is well-known |[3l 

Theorem 4.1. The operator Tq generates a Co-semigroup of isometries (U(t)) t ^o of X given by: 

U(t)f(x, v) = f(x - tv, v), t > 0. 

The operator A = Tq + E with domain 9(A) = 9(Tq) n 9>(L) is the generator of a contractions 
Co-semigroup (V(t)) t ^o given by 

V(t)f(x, v) = exp(-o(v)t)f(x - tv, v), t > 0. 
Let us define now K as the gain operator Q + endowed with the domain of A: 

9{¥) = 9{A), Kf(x,v) = Q + (f)(x,v)=e- 2 f \q ■n\f(x,v i ,)M 1 (w i ,)dwdn. 

Jr 3 xS 2 

It is clear from (fT31) that, for any f € 9(K), 

\ {Af + Kf)dxdv = 0, (4.2) 

Jr 3 xR 3 

while Kf > for any / € 9(K), f > 0. Then, the following generation result is a direct conse- 
quence of [Tl |25l : 

Theorem 4.2. There exists a positive contractions semigroup (Z(f))f^o in X whose generator G is 
an extension of A + K. Moreover, (Z(f))t>o is minimal, i.e. if(T(t)) t> o is a positive Co-semigroup 
generated by an extension of A + K, then T(t) > Z(t)for any t > 0. 
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The natural question is now to determine whether the "formal" mass conservation identity dl.61 ) 
can be made rigorous. Namely, one aims to prove that, for any nonnegative / € X, the following 
holds: 

l|Z(f)/|| = ||/||, W>0. 



The important point to be noticed is the following. If G = A + K, then any function (p € @(G) can 
be approximated by a sequence {(p n ) n c <2){A + K) - Ss(A) such that cp„ — > cp and (A + K)cp n — > 
G<p as n — > oo. In particular, (14.21 ) implies that 

I G(p dxdv = lim I (A + K)(p n dxdv = 0, V<pe 3>(G). 

JjR3xR3 JR3XR3 

Now, for any given initial datum /o € &{G), fo > 0, the solution f(t) - Z(£)/o of (14.11 ) is such that 

4i/wii= r Tj® dxdv = r Gf(t)dxdv=o, 

at JR 3 xR3 JR3 X R3 

i.e. 

11/(011 = ll/oll, Vf>0. 



This means that, if G — A + K, then the solutions to the linear Boltzmann equation (14.11) are 
conservative. On the other hand, if G is a larger extension of A + K than A + K, then there may be 
a loss of particles in the evolution (see O for the matter as well as [2] for examples of transport 
equation for which such a loss of particles occurs because of boundary conditions). Precisely, if 
G i= A + K then there exists /o e X, fo > such that 

HZ(f)/oll< ll/oll for some f > 0. 

This shows that the determination of the domain &{G) of G is of primary importance in the 
study of the Boltzmann equation. This is the main concern of the so-called substochastic pertur- 
bation theory of Co-semigroups O. 

We point out that the question of the honesty of the semigroup governing the Boltzmann equa- 
tion also arises in the study of the space-homogeneous version of the latter equation. Indeed, it is 
the unboundedness of the collision frequency (and consequently that of whole collision operator 
Q) that may give rise to dishonest solutions to the Boltzmann equation. Actually, to prove the 
honesty of the Co-semigroup (Z(t)) t ^o, we will adopt the strategy developed first in [1] and sys- 
tematized in [3]. More precisely, we will show that the gain operator K Minis the assumption of 
El: 

Proposition 4.3. There exists C > such that, for any fixed q > 0, 

ess sup I k(v' ,v)dv' < C. 

\v-Ui\^Q J\v'-Ui\>Q 

Proof. Since the kernel k(v,v') differs from the corresponding one for classical gas, except from 
numerical constants, one can apply mutatis mutandis the technical calculations of H] Section 4.1] 
(see also El p. 329-330]) to get the desired estimate. □ 

As a consequence, one deduces immediately from [1], the main result of this section: 
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Theorem 4.4. The generator G of the minimal semigroup (Z(t))t^o is given by 

G = A + K. 

In particular, the Co-semigroup (Z(t))t^o is honest and 

I Z(t)f(x, v)dxdv = I f(x, v)dxdv for any f e X and any £ > 0. 
Jr 3 xR 3 ' Jr 3 xR 3 
4.2. Consequence on the entropy production. As a direct application of the above result (The- 
orem H31), we give a direct rigorous proof of the linear H-Theorem of Ifl4l . In order to stay in 
the formalism of fl4l . we shall restrict ourselves to the space-homogenous case. Precisely, let Y 
denote the set of functions depending only on the velocity and integrable with respect to velocities: 

equipped with its natural norm || • ||y. For any nonnegative / and g in Y, we define the information 
of / with respect to g by 

Kf\g) = J 3 [f(v)]nf(v) -f{v)\ng{v))dv 

with the conventions OlnO = and xlnO = —oo for any x > 0. This means that the information 
is nothing but the entropy functional H5 for the particular choice of O(s) = s In s. One recalls the 
main result of 11241: 

Theorem 4.5. Let U be a stochastic operator of Y, i.e. U is a positive operator such that WUfWy = 
\\f\\ Y for any f eY,f> 0. Then, 

I(Uf\Ug)<I(f\g) 

for any nonnegative f, g in Y. In particular, if g € Y is a nonnegative fixed point of U then, 

I(Uf\g)<I(f\g), VfeY,f>0. 

According to the results of the previous section, it is not difficult to see that the restriction of 
(Z(t)) t> o to Y is a Co-semigroup of stochastic operators of Y. Since the unique equilibrium state 
M 6 Y is space independent, one sees that (Ti + K)M - and, in particular, 

Zy(t)M =M, \ft> 0. 

Combining this with Theorem 14.5 1 one obtains a rigorous and direct proof of the H-Theorem lfl4l 
Theorem 5.1]: 

Theorem 4.6. Let fo € Y be a given nonnegative (space homogeneous) distribution function with 
unit mass, i.e. ||/olly = 1. Assume that I(fo\M) < oo, then 

j t I(f(t)\M) < 0, (t> 0), 

where f(t) — Zy(t)fo = Z(t)fo is the unique solution to (14.11 ) with /(0) = fo. 

Remark 4.7. Note that, once the conservativity of the solution to the Boltzmann equation as- 
serted by Theorem \4.4\ the above H-Theorem can be proved by usual standard method of kinetic 
theory. However, we insist on the fact that such standard proofs require the solution f(t,v) to be 
conservative and, in some sense, the use of the substochastic semigroup theory. 
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